Numerical simulation of the thermal fragmentation process in fullerene Cqq 
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ABSTRACT 

The processes of defect formation and annealing in fullerene C 60 at T = (4000 - 6000) K 
are studied by the molecular dynamics technique with a tight-binding potential. The cluster 
lifetime until fragmentation due to the loss of a C2 dimer has been calculated as a function of 
temperature. The activation energy and the frequency factor in the Arrhenius equation for 
the fragmentation rate have been found to be E a = (9.2 ± 0.4) eV and A = (8 ± 1) • 10 19 s . 
It is shown that fragmentation can occur after the Cgo cluster loses its spherical shape. This 
fact must be taken into account in theoretical calculations of E a . 
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In spite of intensive experimental and theoretical studies of fullerene C 60 [1], there is still 
controversy concerning some of its physical and chemical characteristics. In particular, a 
significant discrepancy between theory and experiment in the determination of the energy 
D required for the loss of a C2 dimer from fullerene Ceo was repeatedly noted until the late 
1990s. Calculations performed using the most accurate ab initio methods gave the value D 
= (11 - 12) eV [2] (see also one of the most recent theoretical works [3]), which exceeded the 
experimental value of D — (5 - 8) eV by a factor of 1.5 - 2; see references in [4]. Note, however, 
that the fragmentation (dissociation) energy was calculated theoretically by the equation 



where E{C n ) is the energy of the corresponding clusters, while D in the analysis of exper- 
imental data was usually assumed equal to the activation energy E a [4] entering into the 
Arrhenius equation for the fragmentation rate 



where ks is the Boltzmann constant, T is the temperature, and A is the frequency factor in 
s _1 . It should be noted here that the fragmentation energy D and the activation energy E a , 
generally speaking, have a different physical meaning; therefore, their comparison is not quite 
straightforward. 

When E a was determined by Eq. (2), it was commonly assumed that A ~ 10 15 s _1 as in the 
majority of mid-size clusters. More recently, it has been shown in a series of publications [4 - 7] 
that the value of A for fullerene Ceo can be much higher than had been considered previously, 
A ~ (10 19 — 10 20 ) s _1 . Assuming that A = 5- 10 19 s -1 , the authors of [7] reanalyzed the results 
of a great number of experiments and came to the conclusion that E a = (10.5 ±0.2) eV, where 



D = E(C 58 ) + E(C 2 )-E(C m ), 



(1) 



k(T)=A-exp(-E a /k B T), 



(2) 
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the arithmetic mean and the standard error were calculated with "improved" experimental 
values. This new corrected value of E a is in much better agreement with both the theoretical 
value of D [2, 3] and the more recent experimental data [8 - 13]. However, when revising the 
experiments carried out in the past years, the authors of [7], for various reasons, excluded from 
their statistical treatment a number of publications in which the activation energy E a proved 
to be substantially lower than 10 eV even after its recalculation with a different frequency 
factor. Moreover, the choice of the value A = 5 ■ 10 19 s _1 in [7] (or A = 2 • 10 19 s _1 and even 
A = 3.4 • 10 21 s -1 in the subsequent publications [9, 10, 12]) seems to be somewhat artificial. 
This leaves a certain feeling of disappointment and looks like an attempt to fit the experiment 
to theory. 

It is also appropriate to note the following circumstance. In the calculation of the fragmen- 
tation energy D of fullerene C 60 by Eq. (1), it was suggested that both the initial cluster and 
its dissociation product (C 58 cluster) are in their optimal (most energy-favorable) configura- 
tions with the lowest possible energy E{C n ) at a given number of carbon atoms in the cluster 
(n = 60 or 58). This corresponds to so-called adiabatic dissociation. At the same time, under 
real conditions, the C 2 dimer can be detached from a C 60 cluster which is (at least) strongly 
distorted due to an external action and/or thermal vibrations rather than from the perfect 
fullerene. The structure of the C58 cluster remaining after the evaporation of the C2 dimer 
can also differ from the most stable one (standard mass-spectrometry methods can determine 
the number of atoms in a cluster but not its shape). 

From all the above, it follows that a real-time numerical simulation of fullerene C 60 frag- 
mentation is of great interest. It allows the temperature dependence of the cluster lifetime 



to be determined directly and the values of A and E a to be found from Eq. (2) without any 
additional assumptions. The results of simulations of fullerene C 60 dynamics at high temper- 
ature have been presented previously in a number of works. In this case, however, the main 
attention has been paid to studying phase transformations upon heating the cluster [14 - 16] 
or the mechanism of C2 dimer loss from the cluster [17 - 19]. In the present work, all the steps 
of fullerene Ceo fragmentation are studied, starting with the formation of defect isomers and 
ending with the loss of a C2 dimer. We used tight-binding molecular dynamics [20], in which 
the many-body tight-binding potential allowed the electronic subsystem contribution to the 
total energy to be determined more accurately than in the case of simple classical potentials. 
This potential did not require such a large expenditure of computer resources as ab initio 
approaches, thus making it possible to follow the evolution of fullerene C 60 during t « 20 ns 
(while usually t < .1 ns in ab initio calculations). At the same time, for bond lengths, binding 
energies, the energy gap between the highest occupied and lowest unoccupied molecular or- 
bitals (HOMO-LUMO gap), potential barrier heights in the Stone- Wales transformation [21], 
and other characteristics of fullerene Ceo, the tight-binding potential gives values that are in 
good agreement with experimental data and/or ab initio calculations (for more details, see 
[19]). Using Eq. (1) and this potential, we calculated the most energy-favorable configura- 
tions of C 58 and C 60 clusters and obtained D = 11.04 eV, which is also in agreement with 
theoretical data available in the literature [2, 3]. 

The numerical simulation of fullerene C 60 dynamics was performed at a fixed total energy 
of the cluster [22]. This approach to the problem corresponds to a real experimental situa- 
tion in the fragmentation of fullerenes when clusters excited by a laser pulse do not collide 
with each other during the time t < 1 /is required for their fragmentation; that is, they are 
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not in thermal equilibrium with the environment (see [17]). In this case, the cluster tem- 
perature T is a measure of the relative motion of its constituent atoms, that is, a measure 
of the excitation energy after irradiation with a laser. It was determined by the equation 
^k B T(3n — 6) = (E^ n ), where n = 60, and (E^ n ) is the kinetic energy averaged over 10 3 
- 10 4 molecular-dynamics steps (the time of one step is to = 0.272 fs). At the initial time, 
random velocities and displacements were imparted to all the atoms so that the total mo- 
mentum and angular momentum were equal to zero. The forces acting on the atoms were 
calculated by the Hellmann-Feynman theorem using the tight-binding Hamiltonian. After 
that, the classical Newtonian equations of motion were solved numerically. 

Since the total energy, that is, the sum of the kinetic and potential (E pot ) energies, remains 
constant in the process of cluster evolution, its temperature T (as a measure of the kinetic 
energy) decreases when the cluster transforms to a metastable state with a higher value 
of E pot [a lower binding energy E b = 60E pot (l) — E pot (Q0)] and increases in the opposite 
case. The dependence of T on time t is presented in Fig. la at the initial temperature 
Ti n i = (4250 ±15) K. The dependence T(t) exhibits clearly defined steps. They are associated 
with the formation and annealing of defect isomers obtained from fullerene Ceo and from 
each other as a result of the so-called Stone- Wales transformation [23], which represents a 
transposition of two C-C bonds. While carbon atoms in initial fullerene C 60 are located at the 
vertices of twenty hexagons and twelve pentagons isolated from each other, the first Stone- 
Wales transformation at t ~ 2.2 ns gives an isomer containing two pairs of pentagons with 
common sides. Its potential energy is higher than that of fullerene C 60 by AE pot = 1.42 eV 
[19]. As a consequence, the cluster temperature drops stepwise by AT xs 100 K. If each pair 
of neighboring pentagons is considered as one 5/5 defect, then this isomer contains N 5 / 5 = 2 
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such defects (for other defect isomers, the value of E pot is higher; no isomer with N 5 / 5 = 1 
exists for topological reasons). 

At t ~ 2.7 ns, the temperature drops again by AT 100 K. An analysis of the structure 
of the cluster showed that an isomer with N 5 / 5 = 4 forms in this case. Then (at t 3.1 ns), 
one defect and, almost immediately, another one are annealed. As a result, an isomer with 
A5/5 = 2 forms again. After annealing of all defects, the structure of the initial fullerene C 6 o 
is restored at t ~ 4.4 ns. This structure is retained up to t ~ 7.6 ns, after which a sequence 
of transitions starts, leading eventually to the loss of the C 2 dimer at t w 9.5 ns; see Fig. 2a. 
Note that a nonclassical fullerene forms at t 7.9 ns that contains one heptagon along with 
penta- and hexagons, and tetragons appear subsequently. The C 58 cluster that forms after 
fragmentation also represents a metastable isomer of fullerene C 58 . It contains N 5 / 5 = 6 pairs 
of neighboring pentagons, while N 5 / 5 = 3 in fullerene C 58 with the lowest possible energy. 

Though the C 2 dimer in the above example of fullerene C 60 fragmentation is detached 
from a strongly defect C 60 cluster, this cluster still has a closed spherical shell; see Fig. 2a. 
The C 58 cluster forming as a result of fragmentation also has the same shape. However, our 
analysis showed that this is not necessarily the case. We will give an alternative example. The 
dependence T(t) at T ini = (4335 ± 10) K is presented in Fig. lb. The narrow deep minimum 
in T(t) at t 1.6 ns is associated with the formation and rapid annealing (within At ps 0.1 
ns) of an isomer containing a dodecagon along with a heptagon. Actually, the cluster in this 
case opens for a short time, and a large window appears on its surface. The same occurs 
at t ~ 3.6 and at t ~ 4.8 ns. The next time, such a window forms at t ~ 5.2 ns; however, 
now it adjoins two heptagons rather than one. In addition, the cluster also has two pairs of 
neighboring pentagons. The closed surface is not restored up to the loss of the C 2 dimer from 
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the cluster. The C 58 cluster remaining after fragmentation is also opened, see Fig. 2b. 

Since the fragmentation of fullerene C 60 is of stochastic character, it is reasonable to consider 
the mean lifetime r at a certain temperature rather than the fragmentation temperature [24] . 
In order to collect statistical data sufficient for determining the temperature dependence of 
r, we performed calculations at several tens of values of T ini = 4000 — 6000 K corresponding 
to various sets of initial velocities and displacements. We found correlation between T ini and 
the character of the Ceo cluster distortion immediately before fragmentation. The cluster 
before fragmentation is more frequently opened than closed at T ini > 5000 K and the other 
way around at T ini < 5000 K. Although the fragmentation proceeded by means of the loss 
of a C2 dimer in the majority of cases, the evaporation of either one carbon atom, or a C3 
trimer, or a C 4 tetramer occurred several times, as was also observed experimentally [25]. 
With the aim of comparing the results of simulations with experiments on the loss of a C 2 
dimer, we retained for analysis only variants in which fragmentation proceeded through the 
channel C 60 — > C 58 + C 2 . In order to take into account the effects due to the small cluster 
size and the absence of heat exchange with the environment, we used the following equation 
for the fragmentation rate instead of Eq. (2): 



k(T) = A ■ exp 



E a 



kBT(l - Sff) 



(4) 



2CT > 

This equation is derived from the finite-heat-bath theory in the first-order approximation with 
respect to the small parameter E a /2CT [26, 27]. Here, C is the micro canonical heat capacity 
of the cluster. We assumed that C = (3n — Q)ks, where n = 60. Since E a /2CT ps 0.07 at 
T ps 5000 K, the possible difference between C and (3n — 6)/cb does not lead to any significant 
error in the determination of E a from Eq. (4). 

The calculated dependence of the logarithm of the fullerene lifetime r = 1/ A; on is presented 



in Fig. 3, where T* ni = T ini (l — E a /2CT ini ), see Eq. 4. In a rather wide range of T ini = 
(4000 — 6000) K and, correspondingly, r = 1 ps - 20 ns, the results of numerical simulation 
are described well by Eq. (4) with A = const (T). The statistical scatter of the data prevents 
the temperature dependence of A from being reliably determined. However, it may be argued, 
in any case, that this dependence if sufficiently weak in this range of T ini , since the dependence 
of ln(r) on (X^) -1 is approximated well by a straight line, see Fig. 3. The activation energy 
E a and the frequency factor A can be determined by the slope of this line and the point 
of its intersection with the ordinate axis. For the mean values and standard errors of these 
quantities, we obtained E a = (9.2 ± 0.4) eV and A — (8 ± 1) • 10 19 s" 1 (since E a enters both 
into the numerator of the exponent in Eq. (4) and into the renormalized initial temperature, 
it was determined self-consistently by successive iterations). 

Previously, ideas of the loss of a C 2 dimer from a strongly defect C 60 cluster were discussed 
as a hypothesis for the explanation of the discrepancy between the values of D and E a 
[25]. However, these considerations were not confirmed by the corresponding calculations 
(though evidence for the two-step character of fullerene Ceo fragmentation was obtained in 
[28] using the semiempirical Tersoff potential). Note that the value of E a found in our work 
is somewhat lower than the averaged experimental value E a = (10.5 ± 0.2) eV [7]. However, 
if the experiments giving E a < 10 eV were not rejected, as was done by the authors of [7], the 
experimental mean value E a = (8—9) eV would almost coincide with our result. We emphasize 
the fact that the value of the frequency factor A found in our work by direct calculations 
without any assumptions is in reasonable agreement with the values A = (2 — 5) • 10 19 s -1 , 
which have been used in recent years by the majority of authors in the analysis of experimental 
data [7, 9, 11-13], and does not strongly differ from the value A ~ 10 21 s _1 reported in the 
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most recent publication in this field [29]. 

Note in conclusion that it is interesting to calculate the lifetime of fullerene C 60 at lower 
temperatures T = (1000 — 4000) K and possibly to determine the temperature dependence 
of the frequency factor. With molecular dynamics, this task, however, is beyond the possi- 
bilities of the present-day computational technique. For its solution, fundamentally different 
approaches are required based, for example, on Monte Carlo algorithms. 
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Figure captions 

Fig. 1. Time dependence of the fullerene Ceo temperature obtained by numerical molecular 
dynamics simulations. The total energy of the cluster is constant, and the time of one step 
is t = 0.272 fs. The values of T were calculated by averaging the kinetic energy over 10 4 
steps. The initial temperature T ini = (a) (4250 ± 15) K and (b) (4335 ± 10) K. The sequence 
of the formation and annealing of defect C 60 isomers leading to the loss of a C 2 dimer is 
as follows: C 60 -> C 60 (iV 5/5 = 2) -> C 60 (iV 5/5 = 4) -> C 60 (iV 5/5 = 3) -> C 60 (iV 5/5 = 4) -> 
C 60 (iV 5/5 = 2) - C 60 - C 60 (iV 5/5 = 2) - C 60 (iV 5/5 = 4, N 7 = 1) - C 60 (iV 5/5 = 6, N 7 = 1) - 
C 60 (iV 5/5 = 7,7V 7 = 1) -> C 60 (iV 5/5 = 9,iV 7 = 2) -> C 60 (iV 5/5 = 4,iV 7 = l,iV 4 = 1) -> 
C 60 (iV5/5 = 6) -> C 60 (iV 5/5 = 4,iV 7 = 1,7V 4 = 1) _> C 58 + C 2 for (a) and C 60 -> C 60 (iV 5/5 = 
2) -> C 60 -> C 60 (iV 7 = 1, 7V 12 = 1) -> C 60 -> C 60 (iV 7 = l, jV 12 = 1) -> C 60 -> C 60 (iV 5/5 = 2) -> 
C 60 (iV 5/5 = 4,7V 7 = 1) _ C 60 ^ C 60 (iV 5/5 = 4,7V 7 = 1) _ C 60 (iV 5/5 = 2) ^ C 60 (iV 5/5 = 
5,iV 7 = 1) -> C 60 (iV 5/5 = 4,iV 7 = 2,7V 12 = 1) -> C 58 + C 2 for (b), where iV 5/5 , iV 4 , N 7 , and 
A^i 2 are the numbers of pairs of pentagons with common sides, tetragons, heptagons, and 
dodecagons, respectively; see text. 

Fig. 2. Snapshots of a Ceo cluster before and immediately after fragmentation. The time 
interval between two snapshots is At ps 2 ps. For clarity, the distant atoms are not shown. 
The atoms removed upon fragmentation are shown in black. The temperature prehistory of 
fragmentation processes (a) and (b) are given in Figs, la and lb, respectively. 

Fig. 3. Logarithm of the lifetime r (in seconds) of fullerene C 60 before C 60 — > C 58 + C 2 
fragmentation vs. the reciprocal initial temperature (kelvins) including the finite-heat-bath 
correction, see text. The circles are the results of calculations and the solid line is a linear 
least-squares approximation. 
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